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Abstract. The reaction 7p — > pir ^' has been measured with the Crystal Ball / TAPS detectors using the 
energy-tagged photon beam at the electron accelerator facility MAMI-B. Energy and angular differential 
cross sections for the emitted photon 7' and angular differential cross sections for the n° have been deter- 
mined with high statistics in the energy range of the zA + (1232) resonance. Cross sections and the ratio of 
the cross section to the non-radiative process 7P — > pn° are compared to theoretical reaction models, having 
the anomalous magnetic moment k 4 + as free parameter. As the shape of the experimental distributions 
is not reproduced in detail by the model calculations, currently no extraction of k 4+ is feasible. 

PACS. 13. 40. Em Electric and magnetic moments - 13.60.Le Meson production - 14.20.Gk Baryon reso- 
nances with 5 = - 25.20.Lj Photoproduction reactions 



1 Introduction 

The Z\(1232) as a member of the J p = 3/2 + baryon de- 
cuplet acts as the first and only well-isolated resonance 
in elastic pion scattering or pion photproduction experi- 
ments. Its static electromagnetic properties, particularly 
the magnetic dipole moments /i/i , offer important tests for 
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baryon structure calculations in the nonpertubative do- 
main of QCD. Several predictions for the magnetic dipole 
moments of Z\(1232) isobars have been obtained (see ta- 
ble [T]) using constituent quark based models [TJ[2l[3] as well 
as lattice QCD calculations [3]. Recently, also the chiral 
extrapolation of latice results for [i^ including the next- 
to-leading nonanalytic (NLNA) structure of xPT was re ~ 
ported in ref. [5]. Since there are discrepancies between 
the different model calculations, experimental results for 
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Table 1. Predictions of several baryon structure calculations 
for magnetic dipole moments of the Zi(1232) isobars. 



Model 


Mzi++/MJv 








SU(6) 


5.58 


2.79 


0.00 


-2.79 


RQM [1] 


4.76 


2.38 


0.00 


-2.38 


XQSM 


5.40 


2.65 


-0.09 


-2.83 


XBM 


3.59 


0.75 


-2.09 


-4.93 


LQCD [J 


4.48 ± 0.30 


2.24 ± 0.15 


0.00 


-2.24 ±0.15 


NLNA [5] 


4.99 ±0.56 


2.49 ± 0.27 


0.06 


-2.45 ±0.27 
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Fig. 1. Method of investigating electromagnetic moments of 
Z\(1232) baryons using a/1-) Ay 1 transition within the reso- 
nance width and a subsequent decay A — > Nn. 



Z\(1232) magnetic dipole moments are desirable in order 
to constrain this observable at physical quark masses. 

While experimental studies of the N —> A transition 
have provided access to electromagnetic transition mo- 
ments (e.g. magnetic dipole [In^-a and electric quadrupole 
Qn^a transition moments) [6j[7], the static electromag- 
netic properties of the Z\(1232) itself, like ha and Qa, 
are difficult to measure due to its short lifetime of about 
ta — 1CP 23 s. Thus, the experimental method of spin pre- 
cession measurements of fi that have been performed with 
high precision for octet baryons (N, A, S, E) as well as 
for the fi~ decuplet baryon is not possible for short-lived 
states like the 4(1232). 

It has been proposed [5] to determine the magnetic 
dipole moment of the Z\(1232) from the measurement of 
electromagnetic transitions within the resonance width of 
Fa — 120 MeV. In this process the nucleon is excited to 
a A state, which then emits a real photon 7' and subse- 
quently decays into a nucleon and a 7r meson (see fig. [T]), 
hence leading to a Nn^' final state. Because of spin and 
parity conservation, only Ml, E2 and M3 multipoles are 
allowed for this A —¥ A'y' transition. The amplitude for 



Fig. 2. Bremsstrahlung and A-resonant contributions to Nny' 
final states for pion photoproduction (a) and pion scattering 
(b). Only diagrams (a3) and (63) are sensitive to the magnetic 
dipole moments /j,a- 



this process is dominated by magnetic dipole (Ml) radi- 
ation and, therefore, proportional to fj,A, as higher multi- 
pole orders give only very small contributions. The elec- 
tric quadrupole (E2) amplitude vanishes in the limit of 
zero photon energy due to time reversal symmetry [9] and 
furthermore the E2/M1 ratio of about —2.5% from the 
A" — > A transition amplitude indicates that there is only 
a very small quadrupole deformation of the Z\(1232) [B]. 
Magnetic octupole (M3) radiation is suppressed by two 
additional powers of photon momentum with respect to 
the leading Ml order. 

Unfortunately, Nirj' final states can also result from 
bremsstrahlung radiation emitted from incoming and out- 
going protons or charged pions in ttN —¥ Nit or 7 A" — > Nt: 
processes (see fig. [5]). Such contributions are of the same 
order as the A — > Ay' transition and interfere with the 
process sensitive to ha- Thus, any interpretation of exper- 
imental results and determination of the magnetic dipole 
moment will require and rely on an accurate theoreti- 
cal description of all reaction mechanisms contributing to 
Niry' final states. 

This method for the determination of magnetic dipole 
moments of unstable particles was used for the first time 
in the 1970s at the Lawrence Berkeley National Labora- 
tory (LBNL) [TOllTTj and some years later at the Schweiz- 
erisches Institut fur Nuklearforschung / Paul Scherrer In- 
stitut (SIN / PSI) [T2HT3] . Both experiments used inelas- 
tic pion scattering ir + p — > pTT + "f' in order to extract the 
magnetic dipole moment Ha++ of the Z\ ++ (1232) isobar. 
As a result of several theoretical analyses [HllTMTrjlfrT] 
of both data sets, the Particle Data Group [T8] quotes a 
range of fJ>A++ = 3.7 to 7.5 //jv, where /x/v — eh/2m,N is 
the nuclear magneton. These large uncertainties are at- 
tributed to strong contributions mainly from tt + , but also 
from proton bremsstrahlung in this reaction channel and 
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Fig. 3. Model calculations from ref. [9] for energy differential cross sections da/do/ of the emitted photon 7' at different beam 
energies to. Solid lines are predictions for an anomalous magnetic moment of k a + — 3, dashed lines are for n A + = 0. 



model dependencies in the theoretical descriptions of the 
reaction. 

In the case of the Z\ + (1232) a pioneering measure- 
ment of radiative ir° photoproduction 7p — > pir°~f' was 
performed by the TAPS / A2 collaborations at MAMI 
[T9"] in 1999. From this experiment, together with the the- 
oretical description of the reaction from ref. [9], a value 
of 

fi A + = 2.7±^(stat) ± 1.5(syst) ± 3.0(theo) ^ N (1) 

for the magnetic dipole moment of the Zi + (1232) was ex- 
tracted. The experimental precision of this result is lim- 
ited by the poor statistics of around 500 reconstructed 
7P — > pir°j' events in the previous measurement and its 
quite large systematic uncertainties, resulting from an in- 
complete and inhomogeneous angular acceptance. Also, 
the theoretical reaction model from ref. [9] that was used 
to extract the result |T|) for fi A + introduces further un- 
certainties due to model dependencies in the description 
of 7p — > p7r°7'. Thus, a large improvement on the qual- 
ity both of experimental data and theoretical descriptions 
is needed to get a value for fi A + that allows discrimi- 
nation between different baryon structure calculations. In 
order to give a quantitative impression of the experimental 
and theoretical accuracies needed for the determination of 
7 n g- El shows predictions for energy differential cross 
sections dcr/du/ for the emitted photon 7' in the reaction 
jp — > p7r°7'. These are evaluated within the effective La- 
grangian framework of ref. [S] using different values for 
the anomalous magnetic moment k a + . From these curves 
it can be estimated that even in the favourable case at 
oj = 450 MeV and assuming a precise model description of 
contributing reaction mechanisms, a relative experimental 
precision of about 10% is needed to get a determination of 
the magnetic dipole moment fj, A + with an accuracy of ap- 
proximately 0.5 to 0.8 fiN, depending on details regarding 
the sensitivity to [i A + within a particular model frame- 
work. 

This article presents new data on radiative tt° pho- 
toproduction in the Z\ + (1232) energy region, measured 
with the Crystal Ball / TAPS detector set-up and the 
energy-tagged photon beam at the MAMI-B accelerator 
facility in Mainz. This detector system is described in 




Fig. 4. Glasgow photon tagging system at MAMI-B, produc- 
ing a quasi-monochromatic energy-tagged photon beam. 

sect. [5] and offers a much larger acceptance and efficiency 
than the TAPS photon spectrometer used in the previous 
7P — » p7r°7' experiment. The data analysis and identifica- 
tion of 77) — > pir ^' reactions are presented in sect. [3] to- 
gether with a discussion of background contributions from 
7T° and 7r°7r° production and sect. [4] gives an overview of 
different theoretical reaction models for radiative ir° pho- 
toproduction. In sect. [3] the new experimental results for 
various cross sections and other observables are discussed 
and compared to predictions of these model calculations. 
Finally, our conclusions and a discussion of possible future 
developments are given in sect. HJ 



2 Experimental set-up 

The reaction jp — > pir ^/' has been measured at the Mainz 
Microtron (MAMI) electron accelerator facility [3D] using 
the Glasgow tagging spectrometer [2"T1|2"2"] and the Crys- 
tal Ball / TAPS detector set-up. Bremsstrahlung photons 
are produced by scattering the 883 MeV MAMI-B elec- 
tron beam on a 100 /xm thick diamond radiator, while 
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Fig. 5. Crystal Ball / TAPS detector set-up. Additional in- 
ner detector systems (PID, MWPCs) are installed inside the 
Crystal Ball beam tunnel. 

scattered electrons are separated from the main beam 
and momentum-analysed by a magnetic dipole spectro- 
meter and a focal plane detector system made of 353 
half-overlapping plastic scintillators (see fig. 0}. With the 
known beam energy Eq and the energy E e of scattered 
electrons the emitted photon energy ui is given by 

u = E - E e (2) 

The resulting energy-tagged quasi-monochromatic photon 
beam covers an energy range from 208 to 820 MeV at an 
average energy resolution of Aw ~ 2 MeV and a tagged 
photon flux of 2.8 • 10 7 s" 1 . The photon flux is determined 
with an accuracy of about 5% by counting the scattered 
post-bremsstrahlung electrons with the focal plane detec- 
tors of the tagging system and correcting for the loss of 
emitted photons due to collimation. The probability for 
bremsstrahlung photons to reach the target ( "tagging ef- 
ficiency") is periodically measured by a total-absorption 
lead glass counter, which is moved into the photon beam 
line at reduced beam intensity. At 883 MeV electron beam 
energy and with a 3 mm diameter collimator, the tagging 
efficiency is approximately 34%. 

The photon beam impinges on the target cell (length 
4.76 cm) filled with liquid hydrogen and located in the 
centre of the Crystal Ball (CB) detector [231124]. The Crys- 
tal Ball consists of 672 optically isolated Nal(Tl) crystals, 
each read out by an individual photomultiplier tube. Ev- 
ery crystal has the shape of a truncated triangular pyra- 
mid and is about 40.6 cm long, corresponding to 15.7 radi- 
ation lengths. The Crystal Ball covers the full azimuthal 
range and a polar angle range from 20° to 160°, result- 
ing in a solid angle coverage of 93% of 47r. Electromag- 
netic showers are reconstructed with an energy resolution 
of <r E /E = 0.02/( J B/GeV) 1 / 4 and angular resolutions of 
ag = 2° to 3° and er = a e /smd [23]. 

The target is surrounded by the cylindrical Particle 
Identification Detector (PID) made of 24 plastic scintil- 
lation counters aligned parallel to the beam axis. Each 
detector strip is 30 cm long and 2 mm thick. The PID 
measures the differential energy loss of charged particles, 
which together with the total energy deposited in the 



Crystal Ball, is used in a AE/E analysis for separation 
of protons and charged pions. To provide precise track 
reconstruction for charged particles, two cylindrical wire 
chambers (MWPCs) are placed between the PID and the 
Crystal Ball beam tunnel, covering a polar angle range 
from 21° to 159°. Directions of charged particles are recon- 
structed from the intersection points of the particle trajec- 
tories with the detection layers in both chambers, where 
angular resolutions of ag — 1.3° to 2.3° and a^ = 1.4° 
and efficiencies of about 95% for protons and 85% for 
7T^~ mesons are achieved. Angular resolutions for recon- 
structed tracks were determined from the passage of cos- 
mic radiation through both chambers while efficiencies 
were obtained from analyses of 7p — > pn and jp —> nir + 
reactions. 

Polar angles between 4° and 20° are covered by the 
TAPS detector [2"51l2"6] in a configuration with 510 BaF 2 
modules arranged as a hexagonal forward wall at a dis- 
tance of 1.75 m from the target centre. This results in a 
solid angle coverage of the combined Crystal Ball / TAPS 
detector set-up (see fig. [5]) of around 97% of 4n. Each 
hcxagonally shaped BaF 2 crystal has an inner diameter 
of 5.9 cm and a length of 25 cm, corresponding to 12 ra- 
diation lengths. Electromagnetic showers are determined 
with an energy resolution of a E /E = 0.0079/(S/GeV) 1 / 2 + 
0.018 and angular resolutions of less than 1° FWHM [2"6] . 
The two components of BaF 2 scintillation light (fast com- 
ponent with Tf = 0.7 ns, slow component with r s = 620 ns) 
are separately digitised with different ADC gates and can 
be used in a pulse shape analysis (PSA) for particle iden- 
tification. In front of each BaF 2 detector crystal a 5 mm 
thick plastic scintillator tile acts as a veto detector for 
charged particles. 

With this set-up approximately 800 hours of data were 
taken in three run periods during 2004 and 2005. In ad- 
dition about 120 hours of beam time were used for back- 
ground measurements with empty target. Trigger condi- 
tions during all these runs were that the total deposited 
energy in the Crystal Ball is more than approximately 
60 MeV and a combined sector multiplicity in the Crys- 
tal Ball and TAPS is three or more. For the Crystal Ball 
such a sector corresponds to a fixed group of 16 adja- 
cent Nal(Tl) crystals, while TAPS is divided into four 
sectors, each consisting of either 127 or 128 BaF 2 modu- 
les. A sector contributes to the combined multiplicity if 
the deposited energy in at least one of its crystals ex- 
ceeds a threshold of 15 MeV (Crystal Ball) or 25 MeV 
(TAPS). Additionally, an independent trigger condition 
(downscaled by a factor of 48) requiring a multiplicity of 
two or higher was used to record jp — > pir° events for 
consistency checks and calibration purposes. 



3 Data analysis 

3.1 Radiative 7r° photoproduction 

In the beam energy range up to 450 MeV radiative ir° 
production has only a cross section of approximately 70 nb 
which is very small compared to 7p — > pir° and 7p — > 
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Fig. 6. AE/E analysis for charged particles detected in the 
Crystal Ball, using energy depositions Spid and -Etot in PID 
strips and correlated Nal(Tl) clusters. Particles in the upper 
band are accepted as protons, particles in the lower band are 
accepted as -k^ mesons. 



p7r°7r° reaction channels with cross sections of between 
60 to 300 fib for single ir° and up to 1 fib for double ir° 
production. Under some circumstances these reactions can 
fake a pTT 0r y' final state in the detectors due to either loss of 
a decay photon in the double ir° case or false detector hits 
in the single ir° case. In order to separate the 7p — > pn 0r y' 
reaction channel from such background contributions, an 
exclusive measurement of the pir°j' final state has been 
performed. With measured 4-momenta of all particles, the 
reaction kinematics are overdetermined, allowing a variety 
of kinematic checks. 

The first step in the jp — > pir°Y analysis is the selec- 
tion of events with three neutral clusters and one cluster 
identified as a proton. For the Crystal Ball, charged par- 
ticles are identified using azimuthal correlations between 
PID hits and clusters in the Nal(Tl) array. Further sepa- 
ration of protons and charged pions is done by comparing 
the energy loss -Epid in the PID strip and the total de- 
posited energy E to t from the correlated Nal(Tl) cluster 
(AE/E technique, see fig. EI)- Precise direction informa- 
tion for protons detected in the Crystal Ball is obtained 
from tracks reconstructed using the MWPC data. Protons 
detected in the TAPS BaF 2 array are identified using the 
veto detector hits for separating charged and neutral parti- 
cles and the pulse-shape information available from short- 
and long-gate ADC data. 

Measured proton energies are corrected for energy loss 
in the target and inner detector systems as well as for 
different shower propagation compared to photons in the 
detector materials. This is done by applying a correc- 
tion function obtained from exclusive measurements of 
7P — > p7r° reactions, where measured proton energies can 
be compared to expected energies calculated from two- 
body kinematics. 

The 7T° meson is identified and reconstructed from its 
decay into two photons, using the invariant mass m 77 as 



Fig. 7. Invariant 77 mass spectrum for the combination giving 
the smallest deviation from the 7r° mass. Data points represent 
experimental results, while the black lines are from MC simu- 
lations of 7p — > pir° (dashed-dotted), 7p — > pir ')' (solid) and 
7P —¥ pTY°-K° (dashed) , respectively. For comparison with ex- 
perimental data the sum of MC distributions is shown (grey 
solid line) . Vertical lines indicate the accepted range for recon- 
structed 7T° masses. 



the selection criterion. As the reaction jp — > pir°j' leads 
to a three-photon final state, three different permutations 
from two out of these three photons are possible (7172, 
7i73: 7273)- Therefore, photons from the combination re- 
sulting in an invariant mass with the smallest deviation 
from the 7r° mass are assigned to be ir° decay photons, 
with the remaining one being interpreted as 7'. The dis- 
tribution of invariant masses for the best photon combi- 
nation is shown in fig. [7] together with spectra obtained 
from MC simulations of the jp — > p7r°, jp pir°j' and 
7P — » p7r°7T reaction channels. The experimental distribu- 
tion is dominated by background from non-radiative single 
7T° production, where a secondary particle from an elec- 
tromagnetic shower separates and travels some distance in 
the detector before interacting with material and creating 
an additional hit ("split-off") outside the primary cluster. 
In these cases, part of the original photon 4-momentum is 
not taken into account, which explains the shift to smaller 
invariant 77 masses in the case of 7p — > pir° background 
contributions. The low-energy background contribution at 
masses to 77 < 100 MeV originates from electromagnetic 
background processes from the photon beam and, there- 
fore, is not reproduced by simulations of ir° photoproduc- 
tion channels. Such contributions are, however, removed 
by the condition 115 MeV < m 77 < 155 MeV for recon- 
structed 7T° masses. 



Because of the large background contribution from non- 
radiative single 7T° production, candidates for jp — > pir°j' 
events are explicitly tested for the two-body kinematics of 
7P — > pit reactions. For the 7r° meson the missing energy 



E miss {TT°) = E* 



2W 



(3) 
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Fig. 8. Missing energy spectra for reconstructed tv° mesons (left) and protons (right) under the assumption of a 7p — > pw° 
reaction. Data points represent results for jp — } ptt j' candidates; for comparison, the grey solid lines show experimental 
distributions obtained from identified 7p —¥ pir events. Events resulting in missing energies right of the vertical lines are 
interpreted as 7p — > p-K reactions. 
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Fig. 9. Missing mass spectrum calculated from thep7r° system. 
The peak around nij = 0.0182 GeV 2 corresponds to the mass 
of a second 7r°. Data points represent experimental results, 
while the black lines are from MC simulations of 7p — > piv 
(dashed-dotted), jp — > pn ^' (solid) and jp — > pn°iv (dashed), 
respectively. For comparison with experimental data the sum 
of MC distributions is also shown (grey solid line). Vertical 
lines indicate the accepted range for missing 7' masses. 



is calculated, where E* denotes the measured ir° energy in 
the cm. frame, and m p are the tt° and proton masses, 
respectively, and 



W 



1^ + 2wm p 



(4) 



is the total cm. energy at a given photon beam energy uj. 
Values around E m i as {-K°) = MeV indicate events from 
7P — > pir reactions, where the measured 7r° energy cor- 
responds to the expected value calculated from two-body 
kinematics (see fig. [HI)- An analogous calculation is done 
subsequently for the missing energy E m i ss (p) of the proton, 
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Fig. 10. Time coincidence between Crystal Ball / TAPS and 
tagging system. Random coincidences are evaluated from back- 
ground windows Ri and subtracted from events located within 
the prompt coincidence window P. 



and all events not fulfilling both the conditions -E m i ss (7r°) < 
—20 MeV and E m i SS (p) < —4 MeV are interpreted as 
7P — > pir background. 

Double 7T° production events can lead to three-photon 
final states if a decay photon escapes detection either due 
to the incomplete angular acceptance or, in case of a highly 
asymmetric 7r° decay, the photon energy falls below the de- 
tection threshold of about 25 MeV. Such events are iden- 
tified using a missing mass analysis of the detected pir 
system without taking into account the additional pho- 
ton 7'. The mass mx of the third final-state particle is 
calculated from 



m x(P 7r °) = + ni p - Ek — E p ) 2 — (k 



<?tt - P P ) 



(5) 



where w, k denote energy and momentum of the beam 
photon, while the energies and momenta of proton and 
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Fig. 11. Analysis of 7p — > pir reactions. Left: Invariant 77 mass spectrum for events with two photons. Right: Missing mass 
spectrum calculated from the reconstructed it . Solid lines represent experimental results, while the grey shaded distributions 
are from a MC simulation of jp — } piv normalised to the experimental peak maximum. Vertical lines indicate the ranges of ty° 
and proton mass values accepted for further analysis. 



Table 2. Background contributions from simulated single and 
double 7T° production reactions to the number of reconstructed 
events fulfilling the analysis conditions for jp — >• pn ^' . 



Beam 


Total 


Single 7T° 


Double 7T° 


energy to 


events 


background 


background 


350 MeV 


5880 


355 


3 


400 MeV 


11996 


345 


66 


450 MeV 


11839 


224 


1139 



7T° in the final state are given by E p , p p and E„, q^- For 
double 7T° production the mass of the second ir° is re- 
produced, leading to a peak around m| = 0.0182 GeV 2 
(see fig. [S]), while for -fp —> pTr ^' the expected value for 
the squared missing mass is around m\ — GeV 2 , corre- 
sponding to the mass of the photon 7'. Thus, events with a 
squared missing mass in the range —0.014 GeV 2 < m\ < 
0.006 GeV 2 are accepted as jp — > pn°y reactions. 

Finally, remaining random time coincidences between 
the tagging system and the Crystal Ball / TAPS detector 
setup are subtracted. These are determined from back- 
ground events outside the prompt coincidence peak (see 
fig'- HOP ■ Choosing the same width for prompt and random 
windows makes sure that timing conditions are the same 
for all events. Random coincidences are evaluated from 
different time regions Ri to Rio to get a larger sample of 
background events and reduce the statistical errors. 

Several Geant 3.21 based MC simulations were per- 
formed to determine detector acceptance and efficiency 
as well as the remaining background contributions. Kine- 
matic distributions for jp — > pir°j' were generated using 
predictions from ref. [5] for five-fold differential cross sec- 
tions d 5 u/dw'di7 7 'dJ7 T at different beam energies. Simu- 
lations using different values of ka + were performed to 
ensure that efficiencies do not depend on the anomalous 



magnetic moment used as input for the reaction model. 
With these simulations the overall detection and recon- 
struction efficiency has been determined to be about 13.5% 
for an exclusive measurement of the pir°j' final state. The 
simulation of jp — » pir background reactions is based on 
MAID 27 , 28J calculations for differential cross sections 
da/dn^ at different beam energies, while for jp — > pir°ir° 
a phase-space distribution together with a beam energy 
dependence according to the total cross section ov^ from 
ref. [33] is used. These simulations show that the remain- 
ing background contributions are approximately 3% from 
non-radiative single ir° and 4% from double 7r° produc- 
tion. Background from jp p7r°7r° appears mainly at 
higher beam energies around 450 MeV, where it rises up 
to about 10%, while at lower beam energies only a negli- 
gible contribution is observed (see table [2J . The absolute 
background contribution is derived from kinematic distri- 
butions like fig. [5] and is checked for consistency with the 
acceptances for jp — > pir and 7p — > p7r°7r° in combination 
with the particular total cross sections for both reactions. 



For the determination of cross sections the simulated 
background is subtracted from the experimental results, 
leading to a total number of about 27600 reconstructed 
radiative ir° photoproduction events. Cross sections are 
calculated from this number of reconstructed jp — > pir *}' 
events, corrected for detection and analysis efficiency and 
normalised to the number of target protons per cm 2 and 
the incoming photon beam flux. In order to account for 
contributions from the target cell material (125 /im Kap- 
ton) that have been found to be of the order of 2%, cross 
sections are evaluated in the same way for data taken with 
an empty target and subsequently subtracted from the full 
target results. 
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Fig. 12. Invariant masses of 77 pairs for double ir production. 
Only the best combination with minimal \ 2 according to eq. 
((7J) is shown. 



3.2 Single and double tt° photoproduction 
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Fig. 13. Missing mass spectrum for double n production. 
The solid line represents experimental results, while the grey 
shaded distribution is from a MC simulation of jp — ► p7r°7r° 
normalised to the experimental peak maximum. Vertical lines 
indicate the accepted ranges for missing proton masses. 



In addition to the 7^ — > pn 0r y' analysis cross sections for 
non-radiative single and double 7r° production reactions 
have also been extracted from the experimental data. As 
these processes form the major background contributions 
for 7p — > p7r°7', a quantitative understanding of jp — > pn 
and 7p — > pir°ir° is desirable and can be used to check the 
consistency of calibrations and data analysis. 

The analysis of single ir° production has been limited 
to events fulfilling the downscaled trigger condition that 
the sector multiplicity is two or more. From this dataset 
events with two neutral clusters are selected and the ir° 
is identified and reconstructed using the invariant mass 
m 77 of the photon pair (see left panel of fig. ITT|) . The 
identification of the jp —> pir° reaction channel is then 
done by evaluating the missing mass 

m x (7T ) = yj(u> + m p - E^f - (k - q^f (6) 

(see right panel of fig. ITTj) . where u>, k and E^, q v denote 
energy and momentum of the beam photon and the recon- 
structed 7T°, respectively. To determine the total cross sec- 
tion, random and empty target background is subtracted 
and the MAID-based MC simulation of jp — > pn°, as de- 
scribed above, is used to correct the number of recon- 
structed events for detection and analysis efficiency, which 
is about 67% for an inclusive measurement of single 7r° 
production. The total systematic uncertainty for this re- 
action is estimated to be about 6%, which covers uncer- 
tainties of 5% in photon flux determination and 1.8% for 
target length and density as well as 3% for acceptance and 
efficiency corrections. 

In the case of double tt° production events with four 
detected photons are selected and the ir° mesons are re- 
constructed from the best combination of 77 pairs. As 
there are 3 different possible permutations, the combina- 
tion with the best simultaneous reconstruction of both 



pions, i.e. with minimal 

X 2 = ( m ry - "O + ( m $ -m-i^J (7) 

is selected (see fig. IT2j) . Similar to the single ir° case, the 
final identification of the reaction channel jp — > p7r°7r° is 
done using the missing mass (see fig. I13p calculated from 
the reconstructed ir° mesons according to 



mx(7r°7r°) = sj (k + p T - - q£ ] ) (8) 

where k = (uj,k) and pr = 0) denote the 4- momenta 
of the initial beam photon and target proton, while the 
4-momenta of the final state pions are given by — 
(E n ,q n ). After subtraction of empty target and random 
background contributions, the total cross section is cal- 
culated using the detection and reconstruction efficiency 
of about 42% for the inclusive case, determined with the 
phase-space MC simulation described above. Kinematic 
discrepancies between simulation and experimental data 
are taken into account in the systematic uncertainty of 
the efficiency correction, which is estimated to be 5%. To- 
gether with the absolute normalisation uncertainties from 
photon flux and target density, which are the same as for 
the single ir° case, this gives a total systematic uncertainty 
of about 7% for the double ir° production cross sections. 



4 Reaction models 

In order to determine the Zi + (1232) magnetic dipole mo- 
ment from experimental results for ^p — > pir ^ 1 an ac- 
curate theoretical description of all contributing processes 
is required. While first theoretical calculations [3"0"ll3"Tll3"2] 
considered only the Z\-resonant mechanism and, therefore, 
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Fig. 14. Total cross sections for background reactions jp — > pit (left) and 7p — > ptv°-k° (right). Black points represent Crystal 
Ball / TAPS results, black lines show MAID [271128] (solid) and SAID [38] (dashed) predictions for single tt° production. White 
squares and grey triangles are double 7r° results from refs. [39] and [40], respectively. Error bars include both statistical and 
systematic uncertainties for jp — > pn° and statistical uncertainties only for 7p — > pir°-iT . 



could not reproduce experimental results, an extended de- 
scription including several other contributions has been 
presented in ref. [5]. This model was used for the extrac- 
tion of from the previous TAPS / A2 experiment 
in ref. [19]. Its calculations are done within the context 
of an effective Lagrangian formalism, where in the first 
step the Zi-resonant mechanism as well as a background 
of non-resonant contributions (Born terms, vector meson 
exchange) are taken into account for a tree-level descrip- 
tion of the 7p — > pir° process. Then, the additional photon 
j' is coupled in a gauge-invariant way to charged parti- 
cles, leading to the jp — > pir ^' reaction and introducing 
the anomalous magnetic moment as a new parame- 
ter to the jAA vertex, where fi^+ and are related 
according to 

fl A + = (1 + KA+) = {1 + K A +) UN (9) 

with to jv and m A the nucleon and Z\ + (1232) masses, re- 
spectively. 

Since the earlier determination of fj, A + several improve- 
ments have been made to the theoretical calculations. The 
effective Lagrangian approach from ref. [5] has been ex- 
tended in ref. |33] to a unitarised dynamical model of 
jp — > Nttj' reactions (in the following referred to as "uni- 
tary model"), including irN rescattering effects in an on- 
shell approximation for intermediate particles (K matrix 
approximation). Furthermore these rescattering contribu- 
tions are treated in the soft-photon limit ui 1 — » for the 
final 7', where the T matrix for jp — > Nttj' reactions is 
directly proportional to the full T matrix for jp —> Nir 
processes. 

In another approach, descriptions of radiative pion pho- 
toproduction 7p — > Nnj' using chiral effective field the- 
ory (xEFT) have been developed [541135] in order to use 
a more consistent and systematic framework compared to 
the "phenomenological" effective Lagrangian models from 
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Fig. 15. Differential cross section da/df2 n for 7p — > pir°. Black 
points and white squares represent Crystal Ball / TAPS re- 
sults for inclusive (ty° only) and exclusive (p and 7r°) mea- 
surements, respectively, while grey triangles are results from 
ref. gl]. Black lines show MAID [271128] (solid) and SAID [38] 
(dashed) predictions. Error bars include statistical uncertain- 
ties only. 



refs. [5] and [33]. The framework of %EFT allows the cal- 
culation of chiral loop corrections in a consistent way fol- 
lowing quantum field theory, where the chiral symmetry of 
the low-energy strong interaction as well as other general 
principles like unitarity and analyticity are included to 
any given order in a systematic expansion over the energy 
scales and hadronic degrees of freedom. The latest calcu- 
lations in ref. [35] are based on a next-to-leading order 
(NLO) chiral expansion with Z\-isobar degrees of freedom 
(using the 5-expansion power counting scheme described 
in refs. [361137] ) and a next-next-to- leading order (NNLO) 
soft-photon expansion with respect to the emitted 7' en- 
ergy, since this is the order at which fj, A + first appears. 
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Fig. 16. Differential cross sections at different ranges for beam energy uj and total cm. energy W, respectively. Black points 
represent Crystal Ball / TAPS results, white squares are results from ref. [19| . Error bars denote statistical errors, grey shaded 
bands show absolute systematic uncertainties. Black lines are theoretical predictions (using n A + — 2.6) of the unitary model 
from ref. [33] (dashed line) and the X^FT calculation from ref. [35] (solid line). Top: Energy differential cross section dcr/do/ 
for the emitted photon 7'. Middle: Angular differential cross section der/di?y for the emitted photon 7'. Bottom: Angular 
differential cross section da/dfi w for the reconstructed 7r° meson. 



5 Results and discussion 

5.1 Single and double ty° photoproduction 

Before discussing our new experimental results for radia- 
tive 7T° photoproduction we also present total cross sec- 
tions for 7p — > pir a and jp p7r°7r° reactions obtained 
with the data analysis as described in sect. [3] Figure [TJ] 
shows the total cross sections for single and double 7r° pro- 
duction in the energy range from the particular thresholds 
up to a beam energy of 800 MeV. For jp — > pn° the new 
Crystal Ball / TAPS results show good agreement with re- 
cent calculations from MAID [271125] and SAID [35J mul- 
tipole analyses that reproduce all recent data for differ- 
ential and total 7r° photoproduction cross sections up to 
800 MeV. Results for double 7r° production are compared 
to two previous measurements from refs. [39J and [ID; with 
the TAPS and DAPHNE detectors, respectively. The over- 
all agreement with both older data sets is rather good. The 
small differences are attributed to the phase-space-based 



efficiency corrections and covered by the systematic uncer- 
tainties. These are estimated to be 10% for the previous 
TAPS data, 4% for the DAPHNE measurement and 7% 
for our new results. 

In addition, fig. [T2] shows the differential cross section 
da/df^K for both inclusive and exclusive measurements of 
7P — > pir°. In the inclusive analysis, as described above, 
only the 7r° has been reconstructed from its decay into two 
photons, while for the exclusive case the recoil proton has 
also been detected and identified. This limits the angu- 
lar acceptance due to detection limitations for the proton 
(e.g. a minimum kinetic energy of T p ~ 50 MeV in the lab 
frame) to approximately 9„ = 70° to 140° and 9„ > 155°, 
where protons are detected in the Crystal Ball and TAPS, 
respectively. The acceptance gap between both detectors 
is caused mainly by additional material of the inner detec- 
tor readout electronics, which blocks the lab frame range 
from about 9 = 10° to 20° for charged particles. As this 
gives an inhomogeneous and complicated acceptance, any 
remaining charged particles within this angular range are 



S. Schumann el al: Radiative 7r° photoproduction on protons in the Zl + (1232) region 



11 



discarded in the analysis. Both the inclusive and the exclu- 
sive results agree rather well with each other and also with 
a previous TAPS / A2 measurement of jp — > pir from ref. 
[31] except for n° polar angles around 9^ = 130°. These 
deviations are attributed to inhomogeneities in the proton 
acceptance resulting from the complicated geometry of the 
Crystal Ball beam tunnel region. The total discrepancy 
between inclusive and exclusive measurements is about 
4%, which is taken into account as an additional system- 
atic uncertainty for the proton reconstruction efficiency. 
In general, the good agreement of our results with model 
predictions and previous experiments for single and dou- 
ble 7T° production indicates that a reasonable understand- 
ing of the main background contributions to jp —> pn ^' 
reactions has been achieved. 



5.2 Radiative ir° photoproduction 

For radiative ir° production in the jp — > pir°j' reaction, 
fig. [TBI shows energy and angular differential cross sections 
for the emitted photon 7' as well as angular differential 
cross sections for the ir° meson. These cross sections have 
been determined for three different beam energy bins, cor- 
responding to cm. energies W = ^/s starting close to the 
Z\ + (1232) resonance position and increasing to 100 MeV 
above. All cross sections that are not differential in 7' 
energy are integrated over uj' > 30 MeV. Systematic un- 
certainties that are shown as shaded error bands in fig. 
[TB] include uncertainties in absolute normalisation, namely 
5% for photon flux and 1.8% for target density, as well as 
uncertainties in efficiency correction. This covers discrep- 
ancies between simulated and experimental kinematic dis- 
tributions, giving an energy- and angular-dependent con- 
tribution to the acceptance uncertainty which is, on aver- 
age, around 7%. Furthermore, a 4% uncertainty for proton 
reconstruction in exclusive measurements is taken into ac- 
count, as described above. Finally, for subtraction of sin- 
gle and double tt° background a systematic uncertainty 
of 7.5% in the number of subtracted events is assumed. 
Quadratic addition of all contributions then results in a 
mean value of aproximately 10% for the systematic uncer- 
tainties on the measured cross sections for jp — > pir°j'. 

The energy differential cross sections da/dco' in the 
upper panels of fig. \W\ are dominated by a 1/u/ behaviour 
as expected from bremsstrahlung contributions of the ini- 
tial and final protons, although, especially in the highest 
beam energy interval, this shape is slightly modified by 
a broad distribution of emitted photon energies around 
uj' = 80 MeV, caused by interference between proton 
bremsstrahlung and the Z\-resonant process. Both theoret- 
ical models seem to overestimate this effect and predict a 
pronounced peak in this energy range, which is not clearly 
visible in the experimental results. Compared to the ex- 
isting data from ref. the agreement is good, except for 
the low-energy region around emitted photon energies of 
uj' = 40 MeV, where the previous measurement resulted 
in significantly smaller cross sections. 

The angular differential cross sections da/dQ 1 ' show a 
broad distribution with a maximum around photon polar 




to / MeV 

Fig. 17. Total cross section for — >• pn in the energy region 
of the A + (1232). Data points represent Crystal Ball / TAPS 
results, error bars include both statistical and systematic un- 
certainties. Dashed and solid lines show predictions from the 
unitarised effective Lagrangian model |33| and the xEFT cal- 
culation [35], respectively. 

angles of about 110° (sec middle panels of fig. HM . This 
is caused by interference between bremsstrahlung from 
the initial and final protons and the Z\-resonant mecha- 
nism, while a pure Z\-resonant process would produce an 
angular distribution peaking at 90°. Both model calcula- 
tions in general reproduce the shape of the angular dis- 
tributions but overestimate the absolute cross sections at 
higher beam energies, while there is reasonable agreement 
between the old TAPS / A2 measurement and our new 
data. 

Angular distributions dcr/di?^ for the tt° meson have 
been determined for polar angles in the range from 8^ = 
60° to 140° and 6 n > 150° (see lower panels of fig. Hty . 
As in the analysis of 7p — > jot 7' reactions any charged 
particles within the laboratory polar angle range between 
10° and 20° are discarded due to the complicated accep- 
tance in this region. This is comparable to the single 7r° 
case (see fig. [T5]) . as the kinematics for ^p — > pir ^' are 
very similar to jp —> pir° especially for soft photons 7' re- 
sulting from the dominating proton bremsstrahlung. Also 
for this observable, the cross section shape is reproduced 
quite well by theoretical calculations, but the the discrep- 
ancies in the absolute values increase more and more with 
rising beam energy. 

These discrepancies in the cross section predictions are 
mainly attributed to inaccuracies in the description of the 
non-radiative jp — > pit reaction (see fig. I17[) . which is 
connected to the jp —> pir°j' process by the low-energy 
theorem for uj' — > 0. This is particularly important for 
the xEFT calculation from ref. [35] , that does not include 
any vector meson (p, uj) exchange mechanisms. These pro- 
cesses, however, dominate the high-energy behaviour of 
7P — > pTT°, where the A excitation alone yields cross sec- 
tions that are too large and are subsequently reduced by 
the addition of p, uj exchange mechanisms. Some of these 
model dependencies can be avoided by construction of a 
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Fig. 18. Cross section ratio R at different ranges for beam energy u and total cm. energy W, respectively. Black points 
represent Crystal Ball / TAPS results, white squares are results from ref. [19| . Error bars denote statistical errors, grey shaded 
bands show absolute systematic uncertainties. Black lines are theoretical predictions (using n A + — 2.6) of the unitary model 
from ref. [33] (dashed line) and the xEFT calculation from ref. 35 (solid line). 



suitable observable for jp — > pir°j' that does not depend 
on details of the jp — ► pir° reaction description. 

In the soft-photon limit (to' — > 0) radiative 7r° pro- 
duction is completely determined by brcmsstrahlung pro- 
cesses from the initial and final protons and gauge in- 
variance provides a model-independent relation between 
7P — > p7r°7' and non-radiative 7p — » pn° reactions. As 
shown in ref. [12], in this soft-photon limit the three- fold 
differential cross section, after integration over the outgo- 
ing photon angles, is given by 



da/d^Tr 



>o 1 e 2 
to' 2tt z 



da 
dfL 



(10) 



with da/df^Tr being the differential cross section for the 
7jj — > pir° process, and an angular weight function 



W(v) 



1 



2v 



In 



1 



v-1 



1 



with 



4to k 



(k - qf 



(11) 



(12) 



where k, q denote the 4-momenta of the beam photon 
and the 7r° , respectively. Integration over the meson angles 
yields an energy distribution 



da_ 



^° — 



with a weight-averaged total cross section 
= 2 f da 



e 
2^ 



df2 v W(v) 



dn n 



(13) 



(14) 



for the 7p — > pir° reaction. From the low-energy theorem 
of eq. (fT3"j) , as derived in appendix B of ref. [55] , the cross 
section ratio 

R=--'-^j (15) 

cTv dw' 

is defined, with the soft-photon limit value R — > 1 for 
vanishing photon energies lu' — > 0. 



Table 3. Weight-averaged total cross sections ov. Theoretical 
predictions are from the unitary model [33] and the xEFT 
calculation [35]. Experimental values refer to Crystal Ball / 
TAPS results with statistical and systematic errors. 



Beam 


Weif 


;ht-averag 


ed cross section av 


energy ui 


unitary 


X EFT 


experimental 


350 MeV 


56.97 nb 


60.51 nb 


56.18 ±0.06 ±2.26 nb 


400 MeV 


42.17 nb 


51.82 nb 


40.59 ± 0.06 ± 1.64 nb 


450 MeV 


32.31 nb 


49.47 nb 


25.91 ±0.07 ± 1.04 nb 



This ratio is less sensitive not only to experimental 
uncertainties (e.g. photon flux, target density) which con- 
tribute to both da/duj' and ov and hence cancel out in R, 
but also to inaccuracies in theoretical model calculations 
of 77) — > pir° , as this observable depends only on the devi- 
ation from the soft-photon behaviour. Quantitative effects 
of these model dependencies can be estimated through the 
weight-averaged total cross section ov from eq. ([T4"]) for the 
7P — > piT reaction. Table [3] compares predictions for a v 
from both the unitary model [33] and the xEFT calcula- 
tion 35 with experimental Crystal Ball / TAPS results 
obtained from an analysis of jp —> pn at three different 
beam energies u in the considered energy range. While the 
unitarised effective Lagrangian framework |33j gives good 
agreement with experimental results at least in the low 
and medium photon energy range, the ^EFT calculation 
[55] shows discrepancies of about 10% already at low pho- 
ton energies, increasing to nearly a factor of two at the 
highest beam energy. As mentioned above, this is due to 
reaction mechanisms not included in the %EFT expansion 
and, therefore, a well-understood limitation, which does 
not affect the description of Z\-resonant processes. 

With these results for the weight-averaged cross sec- 
tion ov for the non-radiative jp — > pir° reaction and the 
energy differential cross sections da/dui' from fig.[T|)J cross 
section ratios R have been derived according to eq. (IT51) 
for the three different beam energy ranges (see fig. [T8|) . 
Also for this observable, the same value for the anomalous 
magnetic moment results in significant discrepancies 
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between the unitary and xEFT predictions. This shows 
that there is still a remaining model dependence in the 
description of jp — > jot 7', which does not cancel in R 
and, therefore, is not related to details of the 7p — > pn 
reaction model. Such model dependencies, however, limit 
the accuracy of any values for k a + extracted from experi- 
mental results. 

The following considerations and fits to experimental 
results will concentrate on the medium beam energy range 
around oj — 400 MeV, as this is a compromise between 
sensitivity to k a + , which is fairly low at lower beam en- 
ergies, and the agreement between experimental data and 
theoretical calculations, which gets worse with rising beam 
energy. Such a procedure that uses only a limited part of 
our experimental results is not suitable for an unambigu- 
ous determination of the Zi + (1232) magnetic dipole mo- 
ment, but will still give an impression of the statistical 
and systematic precision for p A + that might be achieved 
if improved theoretical models become available. The left 
panel of fig. H9l shows a \ 2 fit of unitary model predictions 
to our new results for R for emitted photon energies up to 
the kinematic limit at u/ ~ 180 MeV with the anomalous 
magnetic moment k a + as free parameter. This fit gives a 
result of k a + = 2.79±gjg(stat) ± 0.39(syst), but the ob- 
tained fit quality of x 2 /ndf = 2.94 (ndf = 28) indicates 
that, even with the restriction to beam energies around 
oj = 400 MeV, the agreement between experimental data 
points and theoretical model calculations in the unitarised 
effective Lagrangian framework [33J is not sufficient for a 
reliable extraction of the Z\ + (1232) magnetic dipole mo- 
ment. As already mentioned in the discussion of the energy 
differential cross sections da/do/ the model overestimates 
a peak structure for energies around oj' — 60 to 100 MeV 
and does not fully reproduce the shape of the experimental 
distribution. 

In the case of the xEFT calculation (see right panel of 
fig- Ull) it is reasonable to limit the fit to emitted photon 
energies up to oj' = 120 MeV [53], which yields a value 
of k a + = 3.95t^2o(stat) ± 0.85(syst) at x 2 /ndf = 2.72 
(ndf = 16). This restriction to low-energy photons 7' is 
motivated by the expansion schemes used in the %EFT 
framework, which require the initial beam photon to be 
in the order of the zi + (1232) excitation energy, while the 
emitted photon 7' has to be soft. Applicability of this 
model is, therefore, limited to beam energies not too far 
away from the resonance point (which gives another rea- 
son for ruling out the higher beam energy bin at oj = 
450 MeV), while the low-energy expansion in oj' becomes 
problematic for 7' energies near the kinematic limit. How- 
ever, also at low photon energies around oj' — 40 MeV 
there still remains a discrepancy between the experimen- 
tal results and the xEFT calculation. Thus, at the moment 
none of these theoretical descriptions seems to give a pre- 
cise description of the experimental results that would be 
needed for a reliable extraction of the anomalous mag- 
netic moment k a + . Furthermore, the fit results for k a + 
obtained from both models are not compatible with each 
other. 



6 Conclusion and outlook 



A new measurement of radiative ir° photoproduction 7J3 — > 
pir°j' in the Z\ + (1232) energy region has been performed. 
The data were obtained with the Crystal Ball / TAPS 
detectors using an energy-tagged photon beam produced 
at the electron accelerator facility MAMI-B. Compared to 
the pioneering TAPS / A2 experiment in ref. [JJ5], from 
which the Z\ + (1232) magnetic dipole moment p A + was 
extracted for the first time, a considerable improvement 
on statistics by nearly a factor of 60 has been achieved. 
This new experiment yielded differential cross sections for 
both the emitted photon 7' and the 7r° with high res- 
olutions in polar angles 9y and 8 V and emitted photon 
energy oj'. These cross sections show a reasonable agree- 
ment with the previous experimental result, but cannot be 
fully reproduced by different theoretical models [531135] for 
7P — > Nttj' reactions. 

Part of this discrepancy vanishes if the cross section 
ratio R from eq. (Til)]) is evaluated, as this eliminates some 
of the model dependencies, but also with this observable a 
reliable and precise extraction of the magnetic dipole mo- 
ment ha+ seems to be prevented by limitations in the ap- 
plicability of the models and their inability to describe our 
experimental data accurately. Thus, substantial progress 
on the model descriptions is required in order to fully ex- 
ploit the present experimental results, as the large im- 
provement in the quality of this data compared to the 
older measurement of radiative ir° photoproduction re- 
veals shortcomings in the models that previously were not 
apparent due to the limited statistics in ref. |19) . In this 
context already different approaches for the on-shell prop- 
agator in the A — > Aj' process, as they have been recently 
studied for three models in ref. [44] . turn out to have sig- 
nificant influence on both shape and magnitude of the 
theoretical predictions for 7p — > pir°'y' cross sections. 

In the meantime also several other observables have 
been suggested that show a larger and more direct de- 
pendence on (J, A +, i.e. not relying on interference effects 
between the Z\-resonant process and bremsstrahlung con- 
tributions. For example, ref. |35] proposes to measure he- 
licity asymmetries with a circularly polarised photon 
beam. These helicity asymmetries are nonzero only for 
three-body final states and vanish exactly in the the soft- 
photon limit oj' — > 0, as -fp — > Nit^' processes reduce to 
7P — > Ntt two-body reactions. Furthermore, JC c irc shows a 
linear dependence on fx A + , where it should be noted that 
this is a model-independent feature because of the low- 
energy theorem for "/p ~ > Nttj' . This model- independent 
determination of p A + would require a measurement with 
a high degree of circular photon beam polarisation in the 
energy region around the Zi + (1232) resonance point. How- 
ever, as the predicted magnitude of these helicity asymme- 
tries is in the order of only about 1%, such a measurement 
would require a further major improvement in statistics 
which seems currently out of reach. 
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Fig. 19. Cross section ratio R for beam energies from ui = 375 to 425 MeV. The dashed and solid lines represent the best 
X 2 fit for the unitary model (left) and the x^FT calculation (right) with k a + as free parameter. Grey shaded bands show the 
variations of the model predictions within the statistical fit errors for n A + , vertical lines indicate the range of emitted photon 
energies used for fitting. 
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